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Abstract. Two possible applications of random decoupling are discussed. Whereas so far decoupling meth- 
ods have been considered merely for quantum memories, here it is demonstrated that random decoupling is 
also a convenient tool for stabilizing quantum algorithms. Furthermore, a decoupling scheme is presented 
which involves a random decoupling method compatible with detected-jump error correcting quantum 
codes. With this combined error correcting strategy it is possible to stabilize quantum information against 
both spontaneous decay and static imperfections of a qubit-based quantum information processor in an 
efficient way. 

PACS. 03.67.Pp Quantum error correction and other methods for protection against decoherence - 
03.67.Lx Quantum computation 
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1 Introduction 

Quantum information processing of many-qubit systems 
requires efficient methods for protecting arbitrary logi- 
cal quantum states against unwanted environmental in- 
teractions or uncontrolled inter-qubit couplings. So far, 
two strategies have been developed for achieving this goal, 
namely dynamical decoupling and redundancy-based quan- 
tum error correction. Dynamical decoupling lj is based on 
the main idea of suppressing unwanted influences on the 
quantum dynamics by appropriately applied determinis- 
tic or stochastic external forces. A disadvantage of this 
approach is that typically residual errors are left. An ad- 
vantage of this method of error suppression is that the 
complete physically available Hilbert space can be used 
for logical information processing. In redundancy-based 
error correction [seed, chapter 10] only a subspace of the 
physically available Hilbert space can be used for logical 
information processing. This subspace has to be selected in 
such a way that possible errors can be corrected perfectly 
by appropriate control measurements and recovery oper- 
ations. Thus, this latter type of error correction typically 
requires a sufficiently large number of redundant qubits 
which cannot be used for logical purposes. This represents 
a serious disadvantage if only quantum information pro- 
cessors with small numbers of qubits are available. Thus, 
in view of nowadays technology it is desirable to develop 
error correcting strategies which combine the positive ad- 
vantages of both approaches in order to achieve the high- 
est possible degree of error suppression with the small- 
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est possible number of redundant qubits. Combined meth- 
ods in which environmental interactions are corrected by 
redundancy-based quantum error correction and in which 
uncontrolled inter-qubit couplings within a quantum in- 
formation processor are suppressed by dynamical decou- 
pling may offer interesting perspectives for the stabiliza- 
tion of few-qubit quantum information processors against 
unwanted errors. 

In the following such a combined error correction scheme 
is presented. The main problem which has to be overcome 
in this context is compatibility. One has to ensure that the 
chosen dynamical decoupling method does not leave the 
code space of the redundancy-based quantum error correc- 
tion method. If this compatibility cannot be guaranteed 
uncorrectable errors may arise. Furthermore, in practical 
applications one is also interested in minimizing the resid- 
ual errors originating from errors during recovery oper- 
ations, for example, and from the dynamical decoupling 
method itself. Currently, no universal solution is available 
for this compatibility and optimization problem. 

Here, first results and solutions of these compatibility 
and optimization problems are presented for a particu- 
lar combined error correcting method which is capable of 
correcting a special case of environmentally-induced deco- 
herence, namely spontaneous decay of qubits. Such decay 
processes may arise from spontaneous emission of photons, 
for example. In our subsequent discussion we concentrate 
on the special but important case that spontaneous de- 
cay processes affecting different qubits are statistically in- 
dependent. This implies that in principle not only error 
times but also error positions can be determined by obser- 
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vation of the individual qubits. In the context of photon- 
induced spontaneous decay this case is realized whenever 
the wave lengths of spontaneously emitted photons are 
significantly smaller than the spatial separations between 
adjacent qubits. Under these circumstances detected-jump 
error correcting quantum codes 0, 0| can be used for an 
efficient redundancy-based correction of resulting errors. 
As an additional source of errors we consider uncontrolled 
inter-qubit couplings within a quantum information pro- 
cessor which can be modeled by Heisenberg-type couplings 
with unknown parameters. It has been demonstrated re- 
cently that such errors can be suppressed significantly 
with the help of embedded dynamical decoupling meth- 
ods in which a deterministic dynamical decoupling scheme 
is embedded into a stochastic one |5|. The deterministic 
decoupling scheme leads to a significant error suppression 
for short times whereas the stochastic decoupling method 
guarantees that the residual error-induced exponential de- 
cay of the fidelity of any quantum state is linear-in-timc 
only and not quadratic as in the case of a pure determin- 
istic scheme [a 0] ■ 

This paper is organized as follows. In Sec. 2 basic con- 
cepts of deterministic, stochastic, and embedded dynam- 
ical decoupling methods are summarized briefly. In Sec. 
3 the error suppressing potential of stochastic dynami- 
cal decoupling methods is studied quantitatively with the 
help of the entanglement fidelity. This latter quantity is 
identical to the mean fidelity for large numbers of qubits. 
Whereas previously derived rigorous lower bounds [6j tend 
to overestimate the fidelity decay of a quantum memory 
significantly this mean fidelity yields reliable quantitative 
estimates on the mean accuracy of achievable stabiliza- 
tion properties for both quantum memories and quantum 
computations. In Sec. 4, finally, a new combined error cor- 
recting scheme is introduced in which a compatible em- 
bedded dynamical decoupling scheme is applied within a 
detected-jump correcting quantum code space. This way 
static random imperfections and spontaneous decay pro- 
cesses can be suppressed simultaneously. 



2 Dynamical Decoupling - Basic Concepts 



Analogous to ([1} let us denote the time evolution due to 
H c (t) alone by U c (t), i.e. 



U c (t) = Texp (~i J H c (t)<It/k\ 



(2) 



(T denotes the time ordering operator.) The time evolu- 
tion in the so called toggled frame U(t) = U%(t)U(t) is 
determined by the Schrodinger equation 



H{t) = Ul(t)H Q U c {t). (3) 



In a simplified scenario the time dependence of H c {t) in- 
volves sequences of <5-pulses applied at times U = iAt, i £ 
N. This time evolution gives rise to so-called unitary bang- 
bang pulses of the form <3 , /(i)<7L i _ 1 ) , whereas the i-th pulse 
is the product of elements g which are chosen from a set of 
unitary operators Q according to a certain selection rule 
f(i). For the corresponding time evolution due to H c (t) 
we obtain 



U c (t)=g m , (i-l)At<t<iAt 



(4) 



and the Hamiltonian H(t) of the toggled frame becomes 
piecewise constant, i.e. 

H{t)=g f m Hogf {t) , (t - I) At < t < iAt. (5) 

In the following several selection rules f(i) and their ad- 
vantages and disadvantages will be discussed. They give 
rise to periodic (or cyclic) decoupling, random decoupling, 
and embedded decoupling schemes. 



2.1 Periodic Decoupling 

Let us chose f(i) = (i— 1 mod n c ) + l so that U c (t) reduces 
to a set of n c periodically repeating unitary transforma- 
tions <7/(i) S Q (compare with[T^). This leads to a cyclic 
unitary control transformation, i.e. 



U c (t + T c ) = U c (t), T c = n c At. 



(6) 



Let us consider a quantum system S whose Hilbert space 
Tls is of dimension N. Typically S represents a quantum 
register containing n q qubits so that N — dim TLs — 2 n " . 
The time evolution U[t) of S is due to a time depen- 
dent Hamiltonian H(t) = Hq + H c (t) which is the sum 
of a static term we would like to suppress (this term may 
represent imperfections in a quantum computer, for ex- 
ample), and a time dependent term which represents the 
control operations we are able to apply. Dynamical decou- 
pling tries to achieve a suppression of Hq. 

According to the Schrodinger equation our system evolves 
in time as 



We can use average Hamiltonian theory (AHT) Q to ex- 
press the unitary time evolution operator of the toggled 
frame at integer numbers n of cycles as 



U(nT c ) = exp {-iHnT c /ti) 



(7) 



Thereby, H is given by the Magnus expansion H = + 
H(l) + #(2) +< . . [|. In the limit of fast control sequences, 
i. e. T c -> and M -> oo for a given time T = MT C > 0, 
it is sufficiently accurate to consider lowest order AHT 
yielding 



J-c Jo 



1 



-At 



J^gjHogiAt. (8) 



U(t) =Texp (~iJ Q H 



(T)dr/h 



In order to achieve the main goal of dynamical decoupling, 
^ ' namely a vanishingly small total time evolution in the 
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(a) fag\ fcg\ ffiflnc 5352 fcgl 91 

■< 1 1 1 1 1 1 

tlme T c + 2At T c + At T c = n c At 2At At 



(b) fsfl fif\ Hf\ ?2r\ fi 
< 1 1 1 1 1 

tlme AAt 3At 2At At 



(c) 3352 92g\ gif2f\gi c g 3 gl g 2 g\ gi fi 
< 1 1 1 1 1 1 

tm "' T c + 2At T c + At T c = n c At 2At At 



Fig. 1. Schematic representation of: (a) periodic dynamical decoupling (PDD) g 6 Q, (b) naive random decoupling (NRD) 
r 6 Q, (c) embedded decoupling (EMD) g 6 Q and r 6 Q' . The time evolution between subsequent instantaneously applied 
unitary operations is assumed to be governed by a perturbing Hamiltonian Ho- 



toggled frame, H^°) has to vanish. This latter requirement 
can be fulfilled by selecting a set Q of unitary operations, 
for example, for which the relation 

5> + i?o<? = (9) 

geS 

holds. Such a set Q of size n c = \Q\ we call a decoupling 
selQ for the Hamiltonian Hq. This first order decoupling 
wor ks perfectly only in the limit of fast control sequences. 

I Viola and Knilll Q calculated a lower bound for the 
worst case fidelity 

F W (T = nT c ) = min|(V>|?7(T)|V>)| 2 (10) 

which is valid for arbitrary periodic dynamical decoupling 
(PDC0) sequences. In particular, they showed 

F^ DD > 1 - 0(T 2 T c V) (11) 

with k — H-ffolh denoting the operator-2 norm of the in- 
ternal Hamiltonian Hq. 

A common procedure for obtaining higher than low- 
est order decoupling is to engineer each cycle in such a 
way that is it time symmetric. This can be achieved by 
traversing an original cycle Q after its completion again 
backwards so that the relation H(T' C — t) = H(t) holds 
with T' c = 2T C . In such a case all terms of odd orders 
in the Magnus expansion vanish and one obtains the im- 
proved bound [l(| 

if u DD > 1 - 0(T 2 T c 4 k 6 ). (12) 

Thereby, SDD stands for symmetric dynamic decoupling. 

A disadvantage of periodic (or cyclic) decoupling is 
the quadratic time dependence of the fidelity decay and 
its dependence on the length of the sequence. As a conse- 
quence for long times and/or large decoupling sets Q the 
accumulating errors might increase quickly. 

1 On how to construct such a set see e. g. [j| 

2 The abbreviations for the different decoupling methods are 
adopted from [ljj. 



2.2 Annihilators 

We defined a decoupling set Q for a Hamiltonian Hq as a 
set of unitary operation g £ Q such that 

E-9 tff o.9 = 0. (13) 

We will call a decoupling set for which the above equation 
holds for all traceless Hamiltonians an annihilator. It was 
shown in [Tl| that such a set contains at least N 2 elements. 
An annihilator can always be chosen to be also a so called 
nice error basis, which is a unitary and orthonormal oper- 
ator basis fulfilling some additional convenient properties 
(for a detailed definition see [12] )• Since the elements of a 
nice error basis form an irreducible projective representa- 
tion of an underlying index group Schur's lemma implies 
the important relation [IJJ 

In our subsequent discussion we will take as annihilator 
the nice error basis 

V = {l,X,Y,Z}® n ", (15) 

which consists of tensor products of the Pauli matrices, 

*-(;!)-*-(:i)- y -(!o)-*-G-i)- 

(16) 

However, any other annihilator would lead to the same 
results. 



2.3 Random Decoupling 

Let us choose f(i) =Er {1, 2, . . . , n c }, i.e. elements /(«) 
are selected randomly from the set {1,2,- • • ,n c } giving 
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rise to randomly selected unitary transformations = 
gf{i) (seeQJ)). To calculate the worst case mean fidelity of 
such a naive random decoupling (NRD) scheme we con- 
sider the average over all realizations of the r% which we 
will denote by E, i.e. 

F W (T = nAt) = minEKV-l^f/CTMl 2 . (17) 

W 

For the moment we call this decoupling scheme 'naive' 
since there are more sophisticated decoupling methods 
which will be discussed in the next subsection. It was 
shown in 6] that this mean fidelity obeys the inequality 

^nrd > j _ ( T AtK 2 ). (18) 

An advantage of NRD methods versus periodic meth- 
ods is the linear time dependence of the fidelity decay and 
its independence of the length of the decoupling set Q. 
However, at the same time the error suppression is weaker 
(compare the factor k 2 of (fT8|l with the factors ft 4 of (jlljl 
or n e of (I12|l). As a consequence of this independence it is 
always possible to draw the random unitary transforma- 
tions from an annihilator, such as Q — V . However, choos- 
ing a particularly adjusted decoupler Q for a Hamiltonian 
Hq may offer advantages if one is restricted to particular 
unitary transformations. 

In addition NRD schemes are well suited for stabiliz- 
ing quantum algorithms against perturbations [7[. This 
aspect will be discussed in detail in section [3l In the lat- 
ter reference we introduced the term Pauli-random-error- 
correction (PAREC) for the stabilization of a quantum 
computation using NRD decoupling with Q = V . 

2.4 Embedded Decoupling 

In embedded decoupling (EMD) schemes two different sets 
of decoupling operations are used, namely a set Q which 
achieves (at least) first order decoupling for a perturbing 
Hamiltonian Hq and an annihilator Q 1 . 

Using decoupling operations in a deterministic way 
from a set Q for PDD [or SDD] we obtain an effective time 
evolution which can be described by exp(— iHT c ). Accord- 
ing to © the Magnus expansion yields i?' ' = for a PDD 
scheme [and = = for a SDD scheme]. We can 
suppress H even further by applying a NRD scheme in- 
volving a set of unitary operations Q' . Thereby, random 
unitary operations from this set are applied inbetween any 
two PDD or SDD] cycles (compare with figure QJ;). 

A bound on the worst case fidelity we can obtained by 
simply taking the NRD bound of (fT8| and replacing At 
by T c and k by a bound on the norm of H. Thus one finds 
the results @ 

F™ D > 1 - 0(TT c 3 k 4 ) (19) 

for an embedded decoupling scheme (EMD) based on an 
underlying PDD scheme and [10( 

F^ EMD > 1 - 0{TT^k 6 ) (20) 



for a symmetric embedded decoupling scheme (SEMD) 
based on an underlying SDD scheme. 

Alternatively one may achieve comparably good re- 
sults by using random path decoupling schemes as dis- 
cussed in In these schemes one uses a PDD (or SDD) 
scheme based on a set Q and in each cycle j the selection 
rule fj(i) — 7Tj((i! — 1 mod n c ) + 1) is used. Thereby, ttj is 
a permutation of the set {1, 2, . . . , n c } chosen at random 
and independently for each cycle. 

3 Stabilizing quantum algorithms against 
static imperfections by naive random 
decoupling 

As mentioned in section [2~3l naive random decoupling (NRD) 
can be used to stabilize arbitrary quantum algorithms 
against static imperfections in a rather simple way. In this 
section basic properties of this stabilization are explored 
quantitatively. As a quantitative measure for the accuracy 
of stabilization the expectation value of the entanglement 
fidelity is used which is approximately equal to the average 
fidelity in the case of high dimensional quantum systems. 
Numerical results demonstrating characteristic stabiliza- 
tion properties of the recently proposed PAREC method 
are presented for the special case of the quantum Fourier 
transform. In particular, they are compared with the re- 
cently proposed stabilization method of Prosen et al. [II] . 

3.1 Entanglement fidelity and average fidelity 

The bound (TT5)) was calculated for the worst case fidelity. 
Typically this bound overestimate the actual fidelity decay 
to a high degree. Therefore, for practical purposes it is 
more convenient to consider other fidelity measures, such 
as the average fidelity 

F{£) = J #<V|£(|^|)|V> ( 21 ) 
or the entanglement fidelity 

F e {S) = {4>\{l®£){m4>\)\4>)- (22) 

Thereby, the integration involved in the definition of the 
average fidelity pip has to be performed over the uni- 
form (Haar) measure on the relevant quantum state space 
with the normalization J dip = 1 and £ denotes a trace- 
preserving quantum operation, i.e. a general completely 
positive map. The pure quantum state \<p) occurring in 
the definition of the entanglement fidelity is a maximally 
entangled state between the quantum system under con- 
sideration and an ancilla system of the same dimension 
N . Furthermore, X denotes the identity operation acting 
on the ancilla system. Accordingly, the entanglement fi- 
delity measures the degree to which the entanglement of 
quantum state is preserved by a quantum operation £. Ap- 
parently, it is independent of the choice of the maximally 
entangled state since any two maximally entangled states 
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are related by a unitary acting only on the ancilla. Both 
fidelity measures are not independent but are related by 



F{£) 



NF ^ ) + 1 =F e{ e )+ o( l -^l). (23) 



N + l 



N 



Thus, in the case of high dimensional quantum systems, 
i.e. TV = 2 n i ^> 1, the difference between both measures 
tends to zero. 

In the following we are mainly interested in the entan- 
glement fidelity comparing a unitary operation U and its 
slightly perturbed version Us- Thus the relevant quantum 
operation £ involves a single unitary Kraus operator K 
which is given by K = U' ■ Us. On the basis of (|2"3"|) in 
the case of high dimensional quantum systems the aver- 
age fidelity is approximately given by the entanglement 
fidelity 

1 



F e (UUu s ) 



N 



tr{U^U s } 



(24) 



which is determined by the absolute square of a fidelity 
amplitude 

A=±tr{rfU t }. (25) 

Typically, the evaluation of this fidelity amplitude is much 
simpler than the direct evaluation of the average fidelity 
(|2~lj) . Therefore, in view of its close relationship to the 
average fidelity our subsequent discussion will concentrate 
on the behavior of the entanglement fidelity mainly. 



3.2 Static imperfections and short-time behavior of 
the fidelity 

Typically, the fundamental unitary transformation U con- 
stituting a quantum algorithm can be decomposed into a 
sequence of n g elementary one- and two-qubit quantum 
gates, i.e. 

U = U, lg U 3 -U 2 -U v (26) 

Let us assume in our subsequent discussion that the quan- 
tum algorithm under consideration involves t interactions 
of such a fundamental unitary transformation U. Such 
quantum algorithms a ppe ar in the context of search al- 
gorithms, for example [1 61 ] . Furthermore, let us focus our 
attention on the case of static imperfection in which the 
perturbing influence on such a quantum algorithms arises 
from unknown (but fixed and time independent) Hamil- 
tonian coupling Hq between the qubits constituting the 
quantum information processor. In the idealized situation 
that the controlled elementary quantum gates Uj, j 



with SHj = HoAt/h. In the notation of (|27| the index 
r = 1, ...,t takes into account that the perturbations could 
be different for successive iterations of the algorithm. How- 
ever, in the case of static imperfections considered here 
this possible complication is not relevant. In Appendix 
A it is shown that after t iterations and for sufficiently 
short times T = tn g At according to second order time- 
dependent perturbation theory the fidelity amplitude A is 
given by 

F e (t) = \A(t)\ 2 

= 1 - 1 E m trU/L/W/;..., • u[ JHU k .. A } 



-2E(*-t) E N tx{U ~ n '■ U L-J mU i~ 1 ' 

T=l j,k = l 

U T -U\ fc dff[/ fe ...i} + 0(5ff 3 ). 

(28) 

Thereby, the first term in the sum of (|28|) describes the in- 
fluence of perturbations occurring in the same iteration r 
and the second double sum describes their influence in dif- 
ferent iterations. This particular form of the short-time be- 
ha yior of the fidelity amplitude has been studied in detail 
bv iFrahm et al.l u\ . In particular, these authors demon- 
strated that whenever an ideal unitary transformation of 
a quantum map U can be modeled by a random matrix for 
sufficiently short times after t iterations the corresponding 
decay of the entanglement fidelity is given 



Fe(t) 



| j 2 t 

t c a t c N 



(29) 



where l/t c w n 2 g tr {Hq} At 2 /N and a denotes the relative 
fraction of the chaotic component of the phase space of 
this map. Furthermore, numerical studies indicate that 
the behavior of higher order terms is such that the fidelity 
decay becomes approximately exponential, i. e. 



Fe(t) 



t 



2 t 2 



' " 1>( t c at c N 



(30) 



3.3 Suppression of static imperfections by increasing 
the decay of correlations 

In special cases in which an ideal unitary transformation 
U is not decomposed into elementary gates we may sim- 
plify l|28p by taking n g = 1 thus obtaining the short-time 
fidelity decay 

F e (t) = 1- 



that the controlled elementary quantum gates Uj, j = t—i , 

1, ...,n g are performed instantaneously after time inter- ^ (£~I T I) — tr{U~ T 6H U T SH} + 0(5H 3 ). (31) 



vals of duration At during which unknown inter-qubit 
couplings perturb the quantum algorithms in the r-th in- 
teraction the ideal quantum gate Uj is replaced by the 
perturbed unitary quantum gate 



U s (t) = eM-i$HJ)U 3 



(27) 



N 

r=-(t-l) 

This expression has been studied previously by iProsenl 
[l]| . It indicates that the faster the decay of the correlation 
function 

^ tr{U- T SHU T 5H} 



6 



Daniel Geberth et al.: Combined dynamical decoupling and detected-jump error correction 



the slower the decay of the fidelity. According to an origi- 
nal proposal bv lProsen and Znidaricl 18| this characteris- 
tic feature of the fidelity decay can be exploited for stabi- 
lizing a quantum algorithm against static imperfections. 
This aspect was investigated in detail by these authors for 
the special case of t = 1 . In this case (l28|) reduces to the 
simpler form 



iterations of a unitary quantum map ? Is it possible to 
suppress all off-diagonal elements of the correlation func- 
tion perfectly ? As explained in the next subsection all 
these questions can be addressed and solved in a rather 
straightforward way utilizing NRD decoupling Q. 



j,k=i 



F e = i- — 1 1 • ,>/// , • i i /, m /'.'....! • o(6ir 

(32) 



1- C( 3 ,k)+0(5H 3 

j;k=l 



iProsen and Znidaricl based their error suppression method 
on the idea to rewrite a quantum algorithm U in such a 
way that for the new gate decomposition the sum over the 
off-diagonal elements of the correlation matrix C (j, k) be- 
comes smaller than for the original gate sequence (thereby 
using possibly even a larger number of quantum gates). 
They considered as an example perturbations of the form 
SH = VS with V being represented by a TV-dimensional 
matrix randomly chosen from the Gaussian unitary en- 
semble (GUE). Thus, on average the matrix elements of 



3.4 Stabilization of quantum computations by NRD 
decoupling 

) 

In this section it is demonstrated that the recently pro- 
posed Pauli-random-error-correction (PAREC) method 0] 
can be understood as a stabilization of a quantum com- 
putation using NRD decoupling. This method allows to 
stabilize quantum algorithms against static imperfections 
by converting the typically quadratic time dependences of 
the resulting fidelity decays (129|) into a linear ones. So far 
only numerical evidence has been presented for this re- 
markable property Q. In this section a rigorous formula 
for the time dependence of the average entanglement fi- 
delity of a PAREC-stabilized iterated quantum algorithm 
is presented. It exhibits explicitly the quantitative depen- 
dence of this fidelity decay on all relevant characteristic 
physical parameters. 

In the PAREC method before each elementary quan- 



V fulfill the condition (V jk V lm ) = 6 jm S kl /N. With this tum § ate U <> * = h-,n g of the r-th iteration (r = 1, t) 

of a unitary transformation U a unitary quantum gate 
r (r) jr) t 



kind of imperfections on average the correlation function 
becomes 



(C(j,k))=S 2 



a?2 (33) 



The l/TV 2 -term comes from the fact that according to our 
assumption of traceless perturbing Hamiltonians also our 
matrices V have to be chosen traceless. (In the case of 
a non-traceless perturbations V this restriction can be 
achieved by the replacement V i— ► V—Ttr V/N). It should 
be mentioned that this latter l/iV 2 -term was not taken 
into account in Ref. [l|| so that these authors investigated 



is applied. Thereby, the unitary gates r| (with 

7g = T) are drawn at random from a decoupling set Q . 
In this particular NRD method (compare with section l2~3f 
Q can always chosen to be an annihilator, such as V (com- 
pare with (|15[) ). Simultaneously the changes on the quan- 
tum algorithm due to these random unitary gates have to 
be compensated by replacing each elementary quantum 
gate Ui of the r-th iteration of the original algorithm by 



r^Uirf . Furthermore, after the last quantum 



gate Un] a final unitary gate r^j 1 is applied. As a result 
each iteration of a unitary transformation U is replaced 
by 2n g unitary quantum gates so that after t iterations 
one obtains the result 

U* = U ...UU = 

' n„ \ u n„ •••'3 '2 U 1 '2 'l u l 'l ' n„ )■■■ 



„(*)t 



the quantity S 2 



In order to demonstrate their idea. lProsen and Znidaricl 
considered the Quantum-Fourier- Transformation (QFT) 
as an example. Typically, this unitary transformation U 
is decomposed into n g = ln q (n q + 2)/2j quantum gates 
which involve Hadamard operations and one-qubit phase 
gates (compare with the left-hand side of figure ^ . In- 
stead, IProse^^nT^nidari^ used a different decomposition 
involving n' g — ln q (2n q + 1)/2J quantum gates. In fig- 
ure [3] the correlation matrix (C{j,k)) is depicted for both 
gate decompositions. Compared to the conventional gate 
decomposition (left) the off-diagonal elements of this cor- 
relation matrix are suppressed significantly by this new 
gate decomposition (middle). Diagonal values are always 
constant, i.e. ((C(j,j)) + l/N 2 )/6 2 = L ^ 

Though of interest this proposal of lProsen and Znidaricl 
leaves important questions unanswered. How can such an 
improved gate sequence be found for an arbitrary quan- 
tum algorithm ? How can this be achieved for repeated 



y...r?Vft.^.^)t.^. r W r Wt). 

(^^..4 l M l)t ^2 (l) ■4 l v< l) t.^.^) 



(34) 



with U' 2k 



Uk, U' 2k _ x = rX'rl'W k = 1, ...,n g . A par- 
ticular PAREC implementation of the quantum Fourier 
transform (QFT) is schematically represented in figure 
[2] for the special case of n q — 4 qubits. Definitely, this 
random application of Pauli operations together with the 
associated change of elementary quantum gates does not 
affect any quantum algorithm. 
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Pi Oi Oa ... ... P«, P{ Pa ^3 PJ 




Fig. 2. Quantum circuit of the Quantum-Fourier- Transformation for n q — 4 qubits (left). The first four gates of the same 
circuit involving the PAREC method (right) 
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Fig. 3. In \((C(j, k)) + l/N 2 )/5 2 ~\ for the QFT with n q — 10 qubits using the usual gate decomposition with n g = 60 gates 
(left), the decomposition by Prosen using n' g = 105 gates (middle) and In [(E{C(j,fc)> + l/N 2 )/5 2 ] using the PAREC method 
with n' s — 2n g = 120 gates (right). 



In order to discuss the stabilizing properties of this 
particular NRD methods let us consider the static imper- 
fections of the previous section. Thus, it is assumed that 
each quantum gate Ui, i — l,...,n g is performed instan- 
taneously while between any two of these quantum gates 
there is a perturbing time evolution originating from the 
static imperfections which can be described by the unitary 
time evolution operator exp(— i5H) with 5H = H^At/h. 
As shown in Appendix A for sufficiently small times such 
a static imperfection results in an average decay of the 
entanglement fidelity of the form 

EF e (t) = l-2tn g ^tv{(5H) 2 } 

- 2t J2 Tr E tr{ C/j r ] 5Hr Uj r^SHr} + 0{5 3 ) 

> l-4tn g ^tr{(6H) 2 } + (D(SH 3 ). (35) 

Thereby, the last inequality can be obtained by recall- 
ing that tr{A^_B} constitutes a Hermitian inner product 
for which the Cauchy-Schwarz inequality applies. Eq. (|35[) 
explicitly exhibits the dependence of the entanglement fi- 
delity decay on the number n g of elementary quantum 
gates and the strictly linear dependence on the numbers 
of iterations of the unitary transformation U. 



Several straightforward improvements of the basic re- 
lation (|35|) are possible. If the durations Atj of the static 
imperfections between adjacent elementary quantum gates 
are not equal but depend on these gates, for example, Q35p 
is modified according to 

EF e (i) > l-t-tr{H^}J2(^j+^t hh ) 2 + 0(6 3 ). (36) 

3=1 

Thereby, t-^ denotes the duration of the perturbation af- 
ter the random pulses. 

It is also possible to apply the random gates not before 
each elementary quantum gate but less often. One random 
quantum gate between each iteration of a quantum algo- 
rithm, for example, is already enough to get rid of the 
terms of ([28} quadratic in t. In this case (|35)l is replaced 
by the inequality 

EF e (t) > 1 - t(n g + l) 2 ^ tv{(6H) 2 } (37) 

at the expense that the term linear in t has a coefficient 
quadratic in the number of elementary quantum gates per 
iteration n g . 

In order to determine the decay of the average en- 
tanglement fidelity of a quantum memory stabilized by 
NRD decoupling we use (|35|) and specialize to the case of 
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one iteration t = 1 of a quantum algorithm consisting of 
n„ = n identity gates with a time At/2 between subse- 



quent quantum gates. This way the time between pulses 
becomes equal to At. Denoting the total interaction time 
between the qubits of the quantum memory by T 
one obtains the result 



n g At 



F e (T = nAt) w 1 - An g ^ tv{H 2 }(At/2) 
= l-TAt±tr{H 2 }. 



(38) 



3.5 Static imperfections modeled by the Gaussian 
unitary ensemble 



of the entanglement fidelity becomes 
E(F e ) = 1 - 2n g 5 2 {\ - N- 2 ) 



28* 



— tr Uj 

N 3 



0(S 3 ) 



> l-4n g S 2 {l- N- 2 ) + 0(S 3 ). 



(42) 



Alternatively this expression can also be derived by sub- 
stituting the relevant perturbing Hamiltonian SH into ()35D 
and setting t = 1 in ([35]) . For the special case of a quantum 
Fourier transform (QFT) the resulting values of E(C(j, k)) 
are shown on the right-hand side of figure [31 In this fig- 
ure they are also compared to the corresponding values 
resulting from the improved QFT proposed by Prosen. 



In this section it is explicitly shown that the PAREC 
method is capable of canceling the off diagonal terms of 
the correlation function (C(j, k)) (|3"3")) perfectly. Let us 
consider the stabilization properties of the PAREC method 
with respect to static imperfections which can be charac- 
terized by traceless perturbing Hamiltonians of the form 
SH = H At/h = (V -XtrV/N) /S with V chosen ran- 
domly from the Gaussian unitary ensemble (GUE). These 
perturbations describe physical situations in which in each 
individual realization of a quantum algorithm the inter- 
qubit Hamiltonian perturbing the dynamics of the qubits 
of the quantum information processor is time independent 
but random. To eliminate such GUE-governed static im- 
perfections we have to choose an annihilator, such as V, as 
a random decoupling set. As a result the fidelity averaged 
over all possible random gate reduces to the expression 



E{F e ) = l-J2 E < C (J> fc )> + °( SH3 )- ( 39 ) 

3,k=l 



According to Eq. 



in the evaluation of expectation 



value of the correlation function we have to consider n' g 
2n g quantum gates. In view of the statistical independence 
of subsequent Pauli operations almost all off-diagonal terms 
of the correlation function vanish, i.e. 



fl 



E(C(j,fc)) 
S 2 



_ i 

^trt/,1 2 



l 

N 2 
1 



if j = k, 

if j even and j = k + 1 







if k even and k 
else 



j + 1 



(40) 



4 Combining dynamical decoupling and 
detected-jump error correction 

In this section a combined error correction scheme is in- 
troduced which allows a restricted form of dynamic de- 
coupling within a detected-jump error correcting quan- 
tum code space. The error suppressing properties of this 
combined error correcting scheme are investigated in the 
presence of spontaneous decay processes and coherent but 
unknown Heisenberg-type inter-qubit couplings within a 
quantum information processor. Simple approximate ex- 
pressions are derived for the expected average fidelity de- 
cay of instantaneous decoupling sequences, which are then 
confirmed by numerical simulations. 



4.1 Spontaneous decay 

4.1.1 Fidelity decay due to spontaneous decay 

If the qubits of a quantum information processor couple 
to their environment, spontaneous state changes can occur 
by emitting energy in form of photons, for example. Such 
processes cause state changes from an excited state of a 
qubit, say |1), to an energetically lower lying state, say 
|0). Let us assume for our subsequent discussion that this 
state |0) is the qubit's stable ground state. The dynamics 
of spontaneous decay processes based on photon emission 
can be described by a master equation in Lindblad form 

HI 



Thereby, it has been taken into account that for all uni- 
tary matrices U the relation 



iV 2 

E|tr{lM| 2 = ^£|tr{C/r,}| 2 = l 



(41) 



holds since the average is performed over all unitary ran- 
dom Pauli gates r 6 V which are elements of an orthonor- 
mal unitary error basis. As a result the expectation value 



n q -l 

p = - l -[H s {t), P ] + ([Li,pL\] + [L iP ,L\]). 



(43) 



Here, p describes the density operator of the quantum reg- 
ister involving n q distinguishable qubits. Their unitary dy- 
namics are described by the Hamiltonian H$ (t) . Typically, 
these dynamics originate from externally applied quantum 
gates. The non-unitary aspects of the time evolution of the 
quantum register are described by the Lindblad operators 



Li 



Kfc|0)ii<l| 



k = 1, . 



(44) 
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which describe the influence of spontaneous emission pro- 
cesses on the qubits of the quantum register. In (|43| it is 
assumed that the distinguishable qubits are separated by 
a distance much larger than the wavelength of the sponta- 
neously emitted photons. As a result, to a good degree of 
approximation the qubits couple to statistically indepen- 
dent reservoirs [l9| . If the decay rates of the qubits Kk are 
all equal, it is possible to construct quantum codes capable 
of correcting the spontaneous-emission induced errors. 

However, before addressing such error correcting quan- 
tum codes let us look briefly at the expected fidelity decay 
of a quantum memory under the influence of such spon- 
taneous decay processes. For this purpose let us consider 
the fidelity 

F(t) = (m\p(t)\m), (45) 

with p{t) obeying Eq.f (T43|) ) with the initial condition /o(0) = 
|-0(O))(^(O)|. If a spontaneous decay event affects qubit 
i an arbitrary state of of this qubit is mapped onto its 
ground state |0)j. As a consequence the resulting n g -qubit 
state is expected to have almost no overlap with the n q - 
qubit state before the decay process. Therefore, the fidelity 
is expected to become very small as soon as an decay pro- 
cess has affected any of the n q qubits. If there are n\ 
qubits in excited states initially, the probability that no 
decay event takes place in a time interval of duration t is 
given by 

pit) = exp(— n\ ■ k ■ t) (46) 

provided the spontaneous decay rate k is the same for 
all qubits. Thus, assuming that the fidelity drops to zero 
as soon a spontaneous decay process has taken place the 
average expected fidelity at time t is approximately given 
by Fit) « p(t) [see also|2(|. 

4.1.2 Detected-Jump Correcting Codes 

Spontaneous decay events as described by ()43|) can be 
corrected by detected-jump correcting quantum codes Q. 
Thereby, all modifications of the ideal dynamics which 
do not arise from any photon decay events are corrected 
passively by an encoding of the logical information in an 
appropriate decoherence free subspace (DFS). This DFS 
is spanned by all quantum state involving a constant num- 
ber of excited (physical) qubits. In order to maximize the 
dimension of this DFS one may choose quantum state with 
n q /2 excited qubits. In order to be able to correct also any 
single error due to spontaneous decay of a known qubit 
by unitary recovery operations one has to construct code 
words by bitwise complementary pairing. In the case of 
four physical qubits, for example, such a one-error detected- 
jump correcting quantum code space is spanned by the 
three orthogonal logical states 

|c ) = 4=(|0011) +e^|1100)), (47a) 
| Cl ) = i(|0101) + e^|1010)), (47b) 
|c 2 ) = -UlOUO) + e ^|1001)). (47c) 



Thereby, the value of relative phase tp can be chosen ar- 
bitrarily. A detailed discussion of these error correcting 
quantum codes and their recovery operations can be found 
in Ref . [|J . In particular, it has been shown that provided 
the error position and error time is known the recovery 
operation 

Ri = Xi Y[ CNOTij Hi (48) 

restores any quantum state perfectly. The time t rec needed 
for this recovery operation can be estimated by the time 
needed for the performance of n q — 1 controlled NOT gates. 

In order to simplify computation on the code space it is 
advantageous to use subspaces of one-error detected-jump 
correcting quantum codes which are equipped with a for- 
mal tensor product structure and have been introduced by 
[2lj | . The main idea is to map each logical qubit onto two 
physical qubits, i.e. |0)l — * |01)p and |1) x, — > |10)p, and 
again supplementing these code words with their bitwise 
complements. In order to be able to distinguish between 
logical two-qubit code words, such as |00) and |11)l, for 
example, on has to introduce two additional physical an- 
cilla qubits. Thus, for two logical qubits, for example, the 
resulting four code words are given by 



|oo) L - 


+ 4(|01 01 01)4 


- 1 10 10 10)), 


(49a) 


|01)l - 


+ 4(|01 10 01)4 


- 1 10 01 10)), 


(49b) 


|10)l - 


+ i(|10 01 01)4 


- 101 10 10)), 


(49c) 


|H)l- 


+ 4(|10 10 01)4 


- 101 01 10)). 


(49d) 


This scheme encodes ul logical 
physical qubits. As outlined in 


qubits in np 
ll| and it 


= 2n L + 2 
is possible 



to construct simple universal logical gates on these code 
spaces. With these quantum codes any single spontaneous 
decay process can be corrected perfectly provided the ap- 
propriate recovery operation is performed perfectly. 

Let us look briefly at the expected average fidelity de- 
cay of a quantum memory which is protected by these 
jump codes and whose qubits decay spontaneously. The 
probability that no spontaneous decay takes place during 
a recovery process is given by 

p'iUec) = exp(-n 9 /2 • k ■ t rec ) (50) 

since on average at most n q /2 qubits are excited during 
the recovery operation. Up to a time t one expects an av- 
erage number n e (t) = ^f-nt oi decay events. Thus, at time 
t the fid elity of the quantum memory is approximately 
given by [2(| 

Fit) = p'(t rec ) n °W = exp(- (^) 2 *rec*). (51) 

Comparing (|5ip with the corresponding uncorrected fi- 
delity decay of (|4"6")) (with n\ = n q /2) one notices a sig- 
nificant improvement provided the number of qubits n q 
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is too large, i.e. Kn q t rec /2 <C 1. However, for large num- 
bers of qubits, i.e. nn q t rec /2 > 1, this method does not 
work. This shortcoming can be remedied by encoding the 
physical qubits in blocks so that each block can correct one 
decay event at a time and by using block-entangling gates, 
such as the ones introduced in Ref. [2(1 HH ■ As a result 
error correction fails only if more than one spontaneous 
decay process occurs in the same block during recovery 
and the fidelity decay can be suppressed considerably in 
comparison with (|5ip . 



4.2 Embedding EMD decoupling inside a one-error 
detected-jump correcting quantum code space 

Let us now investigate to which extent additional errors 
originating from uncontrolled inter-qubit couplings can be 
corrected by embedding a suitably chosen decoupling pro- 
cedure into the code space of a one-error detected-jump 
correcting quantum code. 

For this purpose let us consider a case in which in ad- 
dition to spontaneous decay events also static energy de- 
tunings and nearest neighbor interactions affect the phys- 
ical qubits of a quantum memory. Correspondingly, the 
error Hamiltonian describing these inter-qubit couplings 
is given by 



and cross terms involving X and Y Pauli operations up to 
second order average Hamiltonian theory. But it leaves Z- 
and ZZ-type couplings invariant. Thus, neglecting higher 
order terms caused by XX- and FF-type couplings the 
resulting intermediate average Hamiltonian Hq is given by 



H s (t) = KSiZi+ £ nj KtS K s K j+ i 



K=X.Y,Z j=0 



(52) 



H 



with X, Y and Z denoting the Pauli spin operators. Thereby, 
the quantities 5i denote detunings which may arise from 
an external magnetic field, for example, and the param- 
eters Jkj characterize the strengths of nearest neighbor 
couplings. In (|52|) it is assumed that the qubits are ar- 
ranged on a linear chain. 

The detected-jump correcting quantum codes intro- 
duced above rely on a DFS involving a constant num- 
ber of excited qubits. Therefore, one cannot use arbitrary 
decoupling schemes for suppressing the influence of un- 
wanted inter-qubit couplings, since in general they may 
also change the number of excitations in the code words 
and as a result the code space would be left. Within the set 
of Pauli operations only the Z and swap-type operations 
leave these DFSs invariant. So, we can use a combination 
of Z-based flips and a random decoupling based on the 
symmetric group S n of permutations acting on n q qubits 
(swaps) for performing dynamical decoupling inside a one- 
error detected-jump correcting quantum code space. Thus, 
according to the notions of section [2] this constitutes an 
EMD decoupling method, embedding PDD flip-decoupling 
inside NRD swap-decoupling. In particular, flip-decoupling 
applies the flip-operation 



U F = Z <g) X <g> Z i 



I 



(53) 



after time intervals of duration, say r. This has the effect 
of canceling all 2-qubit couplings of XX- and YY-type 



n q — l 

£ 

i=0 



MiZi 



n q — 2 

£ hJz.jZjZj+i. 

3=0 



(54) 



It should be noted that this type of EM D decoupling method 
is applicable only if the number of physical qubits used in 
the code is a multiple of four with a corresponding odd 
number of logical qubits. Otherwise the flip operation (|53|) 
maps some codewords from symmetric (ip = 0) to anti- 
symmetric (ip — 7r) superposition states which causes the 
recovery to fail. After a time At = mr, where m is an 
even integer, swap-decoupling permutes the logical qubits 
by executing up to n q — 1 pairwise swap operations (trans- 
positions). These swap operations are chosen at random 
in such a way that all permutations of the logical qubits 
are equally probable. Keeping track of the applied per- 
mutations, one knows which physical qubit one has to ad- 
dress in order to access the information that was originally 
stored in, say, the first qubit. As a result of this procedure 
decoupling results from the fact that the code words ef- 
fectively see changing coupling constants between their 
qubits. 



4.2.1 Expression for the expected fidelity of flipswap 
decoupling 

To estimate the fidelity decay of a quantum memory per- 
turbed by Hq and protected by the flipswap decoupling 
discussed in Sec l4.2i we calculate the evolution of a pure 
state, say under the influence of the effective toggling- 
frame error Hamiltonian, i.e. 



U(t = NAt) = T J[ cxpi-iH^At/h) 



with 



j=o 



n — r i n r 3 



(55) 



(56) 



and Pj denoting the toggling-frame Hamiltonian and the 
permutations chosen randomly and independently for each 
time interval. Up to second order perturbation theory after 
time t the quantum state is given by 



JV 



\m) = (i-*Y, H o )At / h 



3=1 



N i-1 



- J2 £ B$>H<P [At/hf \ £ H^ 2 [At/hf) \m) 

(57) 



i=2 j=l 



J'=l 
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Comparing this state with the ideal state \ip) yields the mutations results in 
fidelity 



F, 



c 



N 



i=0 
n q —2 



2 



3=0 



i<j 
n q ~3 

— E 



77, 



(ip\H^\i()))\tp)) (At/K) 



q i=0 
n„-2 



N{At/hf({^H^) p \^)~^\H'oW 



2 (p'+-p'-) E J ».* J «J 

J=2 



z,i+l 



J7 



(62) 



(7V 2 -iV)(ZiV^) 2 ((^l(^) p l^> 



l#ol 



c 2 /7//i 2 



with denoting the probability of finding four bits of a 

(58) 

word with n q /2 excited qubits arranged in such a way that 

applying the four-qubit operator ZZZZ to them yields the 

•j--u / \ i ; • ,i u i u -ui eigenvalue ±1. Furthermore, one finds the relation 

with (. . .)p denoting the ensemble average over all possible 6 



permutations, 



^ = h E pi*Pi- 

q ' Pi65„„ 



(59) 



/ / 

p+ -P- 



3 — 4n q + n q 
The last term one has to calculate is given by 



(63) 



In order to find an explicit expression for the expected 
fidelity we need to calculate the ensemble averages of 77q, 



H' 2 and 



Let us start with the calculation of (H' n )p. Since the 
initial state is an encoded quantum state we are only 
interested in the part of the Hamiltonian which has its 
support on the code space. Therefore, it is convenient to 
introduce the projection operator on this code space 77. 
By combinatorial arguments one can show that the aver- 
age of the term involving the detunings 8i is given by 



77 



i=0 



h5i Zi 



77 = 0. 



(60) 



This is due to the fact that each qubit is influenced by 
any detuning equally often and that inside the code space 
exactly half of the qubits are in excited states. The averag- 
ing of the terms involving ZZ-type interactions is slightly 
more complicated. However, on the code space one obtains 
the relation 



n„ — 1 



E KJ*,iZiZi+i) n = — T7 E hJ ^ n 
i=0 ' q i=0 



1 ' i-n 
= : ci77 



(61) 



with the real number c\ depending on the strengths of 
the inter-qubit couplings J z ^. Similarly, one can calculate 
(77g 2 )p. Cross terms of Z- and ZZ-type vanish after en- 
semble averaging. Remaining terms are proportional to X, 
ZZ or ZZZZ . Averaging these latter terms over all per- 



C3, 



(64) 



and depends on the initial state \tp). This initial state 
can always be expressed as a superposition of the sym- 
metric (ip — 0) basis states, i.e. \ip) = X^'^IV'fc) with 
12k \ a k\ 2 — 1- Therefore, detunings do not matter, since 
the application of any Zi operator to any symmetric basis 
state changes it into its antisymmetric pendant and vice 
versa. It is not difficult to to prove the following upper 
and lower bounds for the positive constant C3 



Ci < c 3 < (c 2 



rig— 1 

E 

i=0 



2x2 



h z 5. 



2h 2 



T E^0- (65) 



l<3 



In general, the value of the quantity C3 depends on the 
quantum state. The upper bound corresponds to the case 
0"i — Sii and the lower bound to the case a* = const. 

Inserting these results into the fidelity (|55|) . the qua- 
dratic-in-time term vanishes since n(H' Q )pII oc 77. With 
t = the fidelity due to coherent errors becomes inde- 
pendent of the initial state and is given by 



F c ~ l-tAt{c 2 -c\)/h 2 



(66) 



with At denoting the time between subsequent swap op- 
erations. Thus, a linear- in-time fidelity decay is obtained 
due to the fact that the EMD decoupling suppresses the 
quadratic-in-time decay resulting from inter-qubit cou- 
plings. This is a main results of this section. This cal- 
culated linear can be used as the leading order of an ex- 
ponential approximation yielding the fidelity decay 



F c ~exp(~tAt(c2-cl)/h 2 



(67) 
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Fig. 4. A typical sequence of instantaneously applied flip-type 
(F) and swap-type (S) decoupling operations and of a recov- 
ery operation of duration t rec : During recovery no decoupling 
operations are executed. Here the time delay between swap 
operations At and the time delay between subsequent flip op- 
erations r are related by At = 2t. 

4.2.2 Implementations and Numerical Simulations 

In a practical implementation of the combined scheme 
discussed above one corrects spontaneous decay events 
with priority by taking into account that recovery op- 
erations must not be interrupted by decoupling opera- 
tions. This compatibility requirement guarantees that no 
additional uncorrectable errors occur during the recovery 
steps. A corresponding typical correction sequence is de- 
picted schematically in Fig. |4j This compatibility require- 
ment may lead to varying time intervals between individ- 
ual decoupling operations. However, for the decoupling 
scheme discussed here this does not present any serious 
problems. Assuming that a successful recovery operation 
recovers a faulty state perfectly while a failed recovery 
approximately results in a vanishing fidelity, the fidelity 
decay of such a combined error correcting scheme is given 
by products of the fidelities of Eqs. (j5Tj) and (|67|) . 

In the following we discuss the fidelity decay in the ide- 
alized case of instantaneous decoupling operations. How- 
ever, the recovery operation is assumed to require always 
a finite total time which grows linearly with the number 
of qubits involved (compare with (|48|) ). 

*rec= (|K-l) + 2)to. (68) 

Therefore, secondary spontaneous decays during a recov- 
ery cause errors. Composite gates, such as CNOT gates, 
are assumed to take 3/2 times as long as the basic gate 
time to- 



4.2.3 Instantaneously Applied Decoupling Operations 

Let us assume that the deterministic flip operations and 
stochastic swap operations are performed instantaneously 
after time intervals of duration r and At = jut with an 
even integer m (compare with FigSJ) . Whenever a spon- 
taneous decay event occurs between two flip operations 
the recovery operation is applied and the decoupling se- 
quence is halted for a time interval of duration t rec . Fur- 
thermore, we are primarily interested in cases of frequent 
decouplings, i.e. t rec 3> t, with spontaneous decay repre- 
senting a weak perturbation, i.e. kt< 1. Thus, according 
to Eqs. (|5Tj) and (f6"T|) after time t = NAt the expected 
fidelity Fjd of a quantum memory whose dynamics are 



governed by the master equation (|43|) with Hamiltonian 
(|52|) is approximately given by 

F ID = exp(-(c 2 - c 3 )(tAt/h 2 ) (1 + nt rec n q /2)) x 

cx P (-(^) 2 Wi) (69) 

with n p h = (n q /2)Kt denoting the mean number of spon- 
taneous decay events taking place during the time interval 
t — NAt. Thus, in order to maximize this fidelity it is op- 
timal to minimize the decoupling time interval At. 

4.3 Numerical Results 

Let us assume that the coupling constants Sj and J ZJ of 
the Hamiltonian (|52[) are time independent but unknown 
and that they are distributed randomly and uniformly 
in the interval [—y/3e,y/3e]. In the following we discuss 
the time evolution of the fidelity of a quantum memory 
with three logical qubits in a coherent state perturbed 
by spontaneous decay processes and by inter-qubit cou- 
plings (compare with Eqs. (|43f and (|52|) . The numerical 
simulations were performed with the help of the quantum 
trajectory method [23( ] . Random elements which had to 
be averaged over several trajectories were the realizations 
of random swap operations on the one hand and times 
and positions of spontaneous decay processes on the other 
hand. 

In the case of instantaneous decoupling operations it 
is best to apply them as often as possible as assumed in 
(f69|) . for example. Fig. [5] shows a 3-qubit and an 8-qubit 
encoded quantum memory. The curves drawn correspond 
to the cases of no decoupling and error correction at all 
(red 0), of decoupling only (blue o), of error correction 
by a jump code only (light blue +), and of the com- 
bined scheme discussed here (pink x). The characteris- 
tic quadratic-in-time fidelity decay due to coherent errors 
is apparent from the curves without decoupling (blue o). 
This decay is superimposed by a (fast) linear-in-time de- 
cay due to (un)corrected spontaneous decay. The more fre- 
quently one applies decoupling operations, the closer one 
gets to the spontaneous-decay-induced limit described by 
(ED- 

5 Summary 

Two applications for random decoupling have been pre- 
sented. 

In section 3 it was shown that random decoupling can 
be used to stabilize quantum algorithms against static 
imperfections in an efficient way. In particular, a simple 
expression was derived for the average fidelity decay of 
a quantum memory emerged which complements already 
known rigorous lower bounds for the worst case behaviour. 

In section 4 it has been demonstrated that errors orig- 
inating from uncontrolled couplings to an environment 
and from couplings among qubits of a quantum informa- 
tion processor can be corrected efficiently by combining 
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Time units 



Fig. 5. Time evolution of the average fidelity whose dynamics 
are governed by Eq. (|43p : 3 physical qubits without any error 
correction or decoupling (red 0), 3 physical qubits with in- 
stantaneous dynamical decoupling only with At = 4to (blue 
o), 8 physical qubits encoding 3 logical qubits with combined 
error correction and instantaneous decoupling with At = Ato 
(pink x), 8 physical qubits encoding 3 logical qubits without 
decoupling (light blue +) . (The time is plotted in units of the 
elementary gate time to.) The full lines represent the corre- 
sponding approximate analytical expressions given in the text. 
The parameters are Kto = 10 -3 , eto = 10 -4 . Statistical fluctu- 
ations originate from the sampling of the permutations. In the 
curve for n q = 8, for example, 500 random permutations are 
chosen per quantum trajectory over a total of 100 trajectories 
which is only slightly more than the size of the permutation 
group (8! = 40320). 

redundancy-based quantum error correction with dynam- 
ical decoupling methods. In the particular example pre- 
sented an embedded dynamical decoupling scheme was 
embedded into a one-error correcting detected-jump cor- 
recting quantum code. This way both spontaneous decay 
processes and uncontrolled inter-qubit couplings could be 
suppressed significantly. An important aspect which has 
to be taken into account in the development of such com- 
bined error suppression methods is compatibility. It has 
to be ensured that the dynamical decoupling operations 
used do not leave the code space of the redundancy-based 
quantum error correcting quantum code. Typically this 
compatibility requirement constrains the structure of the 
admissible dynamical decoupling operations. 
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financial support by project MSM 6840770039. Stimulating dis- 
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A Appendix 

In this appendix the perturbative short-time approxima- 
tion of the fidelity amplitude used in (j2"5|) is derived. Let 
us consider a quantum algorithm iterating an ideal uni- 
tary transformation U — U n ■ ■ ■ U3 ■ U2 ■ U\ t times. If this 
ideal time evolution is perturbed by static imperfections 
originating from Hamiltonian inter-qubit couplings the j- 



th quantum gate of the r-th iteration of U is replaced by 
the perturbed unitary quantum gate 

Uj 1 ^ expi-iSHj^Uj exp(-iSHJ r ). (70) 

In Eq. (|28|) . for example, the special case SHJ t = SH = 
Ho At and SHJ r — is considered. The index r in (fTTTj) 
takes into account that perturbations may be different in 
successive iterations of the unitary transformation U. In a 
second order expansion with respect to SH J ; and SH J r after 
t iterations the fidelity amplitude A (|25|) of the iterated 
perturbed quantum algorithm is given by 

' " s 1 1 

= 1 - E E E n ( l tT{SH ^ } 2 tr «^ 2 >) 

t n g j-1 

-EEE^H^ow(fc)} 

t=1 j=2 k=l 

+ tv{5HJ l {j)SHl r {k-l)} 
+ tr{fflj r (j-l)fflj,(fc)} 
+ tv{SHj r (j-l)6Hl(k-l)} 

t n g 

-EE^ tr ^>^^} 

T=lj=l 

t ri-1 n g ^ 

~ E E E 7Y 

Tl=2 T 2 =l j,k=l 

tr{U T2 - Tl 5Hjl(j)U Tl - T2 SH^(k)} 
+ tr{U T2 - Tl SHjl(j)U Tl - T2 SHl 2 r (k - 1)} 
+ tr{U T ^5H^{j - l)U Tl - T2 SH T k 2 (k)} 
+ tT{U T2 - n SHj^(j - l)U Tl - T2 SHH(k - 1)} 
+ 0(SH 3 ) 

(71) 

with 

SHj p (i) = uM ...ut-mr p .Ui... UiUx 

= ul A 5Hj p U i ... 1 . 

The term linear in the perturbing Hamiltonian SH van- 
ishes if all Hamiltonians involved arc tracclcss. Note that 
all the terms of (fTTj) involving tr{-} terms are real val- 
ued so that up to second order the fidelity F e (t) — \A(t)\ 2 
is simply obtained by multiplying all these terms of A(t) 
with a factor of magnitude two. In the special case SHj{ = 
SH = H At, SHJ r = dZU) gives rise to Eq.®- 

Setting SHI = 5HJ r = r^SHr^ Eq. dTTJ) yields the 
time evolution of (|34p . Performing an average E over all 

random unitary operations {t^ } from the decoupling set 
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Q one obtains the expectation value of the amplitude 
EA(t) = l-^^ltr{(rf ) t dffr «)2 } 

r=l j=l 

r=l j=l 

+ C((5 3 ). (73) 
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